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ABSTRACT 

Optical spectra of 25 H II regions in the inner two kpc of the M33 disk have been obtained with the GMOS 
spectrograph at the Gemini North telescope. The oxygen abundance gradient measured from the detection of 
the [O in] A4363 auroral line displays a scatter of approximately 0.06 dex, a much smaller value than recently 
reported by Rosolowsky & Simon in this galaxy. The analysis of the abundances for a large sample of H II 
regions derived from the R23 strong-line indicator confirms that the scatter is small over the full disk of M33, 
consistent with the measuring uncertainties, and comparable to what is observed in other spiral galaxies. No 
evidence is therefore found for significant azimuthal variations in the present-day metallicity of the interstellar 
medium in this galaxy on spatial scales from ~ 100 pc to a few kpc. A considerable fraction of M33 Hll 
regions with auroral line detections show spectral features revealing sources of hard ionizing radiation (such as 
He II emission and large [Ne III], [O III] line fluxes). Since R23 is shown to severely underestimate the oxygen 
abundances in such cases, care must be taken in chemical abundance studies of extragalactic H II regions based 
on this strong-line indicator. 

Subject headings: galaxies: abundances — galaxies: ISM — galaxies: individual (M33) 



1. INTRODUCTION 



The occurrence of radial metallicity gradients in spiral 
galaxies is presently understood within the inside-out forma 
tion paradigm as due to disk growth by mean s of radially 



dependent gas infall and star formation rate ( Matteucci & 
Francois 1989; Boissier & Prantzos 1999). Recent age esti- 
mates of resolved stellar populations as a function of galac- 



populat 

tocentric distance in nearby galaxies has provided convinc- 
ing evidence for the i nside-out evolution o f disk galaxies 
( [Williams eiTaTl|2009| |Gogarten et~aL]|2010| ). At the same 
time, the study of present-day metallicity gra dients in galax 



ies from Hll region spectroscopy (Vila-Costas & Edmunds 
[r992HZantsky et al.|[T994| jKennicutt et al.||2003j ) continues 



to provide important empirical constraints for the complex set 
of pa rameters that describe galacti c chemical evolution mod- 
els flYin et al.|2009l|Fu et alJ2009] ). 

Azimuthal variations in the chemical composition of exter- 
nal spiral galaxies are observationally much less constrained 
than radial gradients. The mixing of the interstellar medium 
(ISM) occurs over a range of timescales, that grow with the 
extent of the spatial scale considered. On scales between 0. 1 
and 1 kpc the effects of cloud collisions and expanding su- 
pershells, powered by stellar winds and supernova explosions, 
and combined with galactic differential rotation, lead to heavy 
element dispersal and mixi ng of the ISM on r elatively short 
timescales, around 10 8 yr ( |Roy & Kunth|l995| ), which helps 
to explain the high degree of chemical homogeneity observed 
in the ISM of galaxies (Scalo & Elmegreen 2004; Edmunds 
2005| ). In the case of the spiral galaxy M101, |Kennicutt &| 



1 Based on observations obtained at the Gemini Observatory, which is 
operated by the Association of Universities tor Research in Astronomy, 
Inc., under a cooperative agreement with the NSF on behalf of the Gem- 
ini partnership: the National Science Foundation (United States), the Sci- 
ence and Technology Facilities Council (United Kingdom), the National 
Research Council (Canada), CONICYT (Chile), the Australian Research 
Council (Australia), Ministerio da Ciencia e Tecnologia (Brazil) and Min- 
isterio de Ciencia, Tecnologia e Innovacion Productiva (Argentina). 



Garnett ( 1996 ) concluded that the local scatter in metallicity is 
significantly smaller than the 0.1-0.15 dex dispersion derived 
from the oxygen abundances of 41 H II regions, a conclusion 
that is in agreement with the expected level of inhomogeneity 
of the ISM from hierarchical models of star formation (~0.05 
dex, |Elmegreen| 1 998]) . 

More recently, |Bresolin et aT] ( |2009| ) measured the oxy- 
gen abundances of 28 Hll regions in NGC 300 using the 
[Om]A4363 auroral line as a temperature diagnostic, and 
found a scatter in t he radial gradient of only 0.0 5 dex. On the 
other hand, in M33 Rosolowsky & Simon ( 2008) found a sub- 
stantial intrinsic fluctuation in the H II region oxygen abun- 
dance, ~ 0. 1 1 dex, in addition to the scatter that is simply due 
to observational errors. In particular, in the central two kpc of 
the M33 disk the dispersion in O/H is 0.2 dex, with a peak-to- 
peak variation of 0.8 dex over a spatial scale of about 1 kpc. 
This result appears difficult to reconcile with the short mixing 
timescale of the ISM and with the small scatter measured in 
the present-day metallicities of other nearby spirals. 

In this paper a new sample of H II regions located in the in- 
ner two kpc of M33 is presented, with the goal of verifying 
the subst antial local dispersion in the o xygen abundance de- 
tected by |Rosol^s^^Sim^|200§ =RS08). Measuring 
accurate Hll region abundances m the central regions of spi- 
ral galaxies is critical for radial gradient studies, but it is gen- 
erally made difficult by the enhanced cooling resulting from 
the increased metal line emission, which hinders the detec- 
tion of the weak electron tempe rature diagnostic lines, such 
as [O in] A4363 (Bresolin 2008). The observations presented 
here alleviate this difficulty, thanks to high-quality detections 
of the [Olll]A4363 line in Hll regions with galactocentric 
distances as small as 0.2 kpc (Sect. 2). The chemical abun- 
dance analysis is carried out complementing the new data 
with existing samples of H II regions, using both the auroral 
line information and the R23 strong-line abundance indicator 
(Sect 3). The systematic effect on R23 -based abundances in- 
duced by the hard ionizing radiation field detected in several 
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of the M33 H II regions is discussed in Sec. 4. 

2. OBSERVATIONS 

Spectra of H II regions in M33 were ob tained with the Ge m- 
ini Multi-Object Spectrograph (GMQS, [Hook et al.|20Q4| ) at 
the Gemini North facility. The targets were selected from 
narrow-band Ha images of two 5.'5x5f5 GMOS fields in 
the inner parts of the galaxy, centered approximately 2. / l 
(0.83 kpc deprojected distance) W and 5 f A (1.67 kpc) S of the 
center. The spectroscopic data were acquired in queue mode 
on October 17, 2009, using two multi-object masks, one per 
field, with 1. "2- wide slits. The seeing conditions were ~0."7 
during the observations of the W field, and ~1"2 for the S 
field. In order to minimize the effects of the differential atmo- 
spheric refraction the data were acquired at airmasses smaller 
than 1.16. Three 1800 s exposures were secured for each of 
the two fields using the B600 grating, which provided spec- 
tra covering the 3500-5100 A wavelength range at a spectral 
resolution of ~5.5 A. For some of the targets, depending on 
their spatial distribution, the spectral coverage extended up to 
-6000 A. 

iRAlQ routines in the gemini/gmos package were used 
for electronic bias subtraction, flat field correction and wave- 
length calibration of the raw data frames. Observations of the 
spectrophotometric standard star Wolf 1346 yielded the flux 
calibration. The final version of the spectra was obtained by 
averaging the three spectra corresponding to each individual 
slit. The H II region sample comprises 25 objects, whose lo- 
cations are shown in Fig. [T] The celestial coordinates of the 
targets are summarized in Table [T] (where objects are listed in 
order of decreasing declination), together with their galacto- 
centric distances in kpc, dep rojected adopting an inclination 
angle for M33 of 56 degrees (Zaritsk y" et al.|1989|), a position 
angle of the major axis of 23 degrees ( |de Vau couleurs et al. 
1991), and the center coordinates measured by Massey ~t al. 
(1996). A distance of 840 kpc from |Freedman et al.| (|2001) 
was assumed. The identification of the subcomponents of the 
Hll regions given in column (5) of Table [T] was taken from 
Hodge etal.| ( [2002l ). 



Table 1 

Observed H II region sample 



2.1. Line intensities and oxygen abundances 

The emission line intensities were measured with the 
splot program in IRAF by integrating the fluxes under the 
line profiles. The line intensities were corrected for interstel- 
lar reddening by assumin g a case B intr insic H7/H/3 ratio of 

( 1979| reddening law. The 
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0.47 at 7; = 10 4 Kandthe_ 

effects of the underlying stellar populations on the strength 
of the B aimer lines were accounted for by requiring that the 
H7/H/3 and HJ/H/3 line ratios provided the same value for the 
extinction. A median equivalent width of 2 A was found for 
the absorption component, and was applied to the whole sam- 
ple. The resulting reddening-corrected emission line fluxes, 
normalized to Hf3 = 100, are presented in Tabled The errors 
account for uncertainties in the flat fielding, theflux calibra- 
tion, the level of the continuum in proximity of each line and 
the logarithmic extinction coefficient c(H/3). 

The electron temperature (T e ) diagnostic line [O in] A4363 
was measured for eight of the 25 Hll regions. The nebular 

2 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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Note. — Units of right ascension are hours, minutes and seconds, and units of dec- 
lination are degrees, arcminutes and arcseconds. Col. (1): Hll region identification. 
Col. (2): Right ascension. Col. (3): Declination. Col. (4): Deprojected galactocentric 
distance in kpc. Col. (5): Identification from Boulesteix et al. 1 1974 = B), Hodge et al. 
(1999] = H), rCourtes et al.| ( 1987 = C) and Magrini et al. 1 2000, = M). O bjects with a 
* symbol in column (5) are in common with Rosolowsky & Simon (2008). 

package in IRAF was used to derive T e for the ++ -emitting 
region, adopting an electron density of 10 2 cm" 3 and the same 
atomic parameters used by |Bresolin et al.|pQ09| ). The temper- 
ature in the + -emitting region was estimated from the rela- 
tion T[0 II] = 0.7 T[0 in] + 3000 K ( |Garnett|1992| |. The ++ 
temperatures and the computed O/H abundance ratios are pre- 
sented in columns (10) and (11) of Table [2] respectively. 

Ten objects in Table [j] are in common with the work by 
RS08 (these are identified by the * symbol in column 5). 
Moreover, NGC 595 (object 1 in Table [T]) has also been re 
cently observ ed with high resolution spectroscopy by Esteban 
et al.| ( |2009| ). A comparison with the line fluxes contained 
in these two publications is presented in Fig. [2] It can be 
seen that the agreement is generally acceptable, except for tar- 
get 15 (C lAb), for which RS08 measured a much stronger 
[Om]A5007 flux than the value reported here. It appears 
likely that the slits used by RS08 and in this work covered 
different portions of this high-excitation nebula. Fig. |2j: in- 
dicates that R23, which can be used as a metallicity indica- 
io^\ is generally consistent between the different works. The 
[O III] A43 63 -derived oxygen abundances in Fig. [2]l can dif- 
fer between this work and RS08 by significant amounts, up 
to 0.4 dex. The logarithmic extinction c(H/3) was found to be 
systematically higher by an average of 0. 19 dex in comparison 
with the RS08 values, for objects in common. The compari- 
son with the fluxes and the oxygen abundance published for 
NGC 595 by |Esteban et^lpO^l full symbols in Fig.[2]) is 
excellent. The c(H/0 value determined for this target is also 



5 R 23 = ([0 11] A3727 + [O ill] AA4959, 5007) / H/3 jPagel etal.|l97? 
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Figure 1. Location of the H II regions studied in this work on a narrow-band Ha: image of M33 (taken from the Local Group Galaxies Survey, Massey ~et al.| 
2006 ). Targets where the [O in] A4363 line was detected are marked with triangles. The galaxy center is shown by the diamond symbol. The objects B 72 and 
BWwere observed with the Subaru telescope by Bresolin et al. (2010). 



in good agreement with their work (0.46 vs. 0.39). 

3. THE RADIAL ABUNDANCE GRADIENT AND ITS SCATTER 

The radial abundance gradient in M33 has been investi- 
gated from H II region spect roscopy by several authors in the 
past f ew years (see RS08, Magr ini et al.|2010| , Bresolin et al. 
2010 and references therein). The large number of available 
H II regions with published line flux data is utilized here to in- 
fer the scatter in abundance at a given galactocentric distance. 
The new Gemini sample was designed to increase the number 
of observations near the galactic center, where the detection of 
the auroral lines is most critical for the radial abundance gra- 
dient, although problematic even with 8m-class telescopes. 

Line fluxes for 83 H II regions with [O III] A4363 detections 
(in addition to [O II] A3727 and [O HI] A5007) have been col - 
lected from the following papers: [K witter & Aller (1981), 
Vflchez et al.| ( [T988]), |Crockett et al.| O006| ), RS08, |JVfagrini 
et al.| (2010) and Bresol met al.| ( |2010| Subaru observations). 
The oxygen abundances were recomputed from the published 
reddening-corrected emission line fluxes with the same proce- 
dure and atomic data used for the Gemini H II region sample. 
The resulting radial abundance gradient is shown in Fig. |3j 
where different symbols are used for different data sources. 
The large scatter of the data points in the diagram, already 



noticed by RS08 for their sample, is evident: the rms scatter 
around the least square fit (dashed line) for the full sample is 
0.17 dex, apparently larger than what can be accounted for by 
the measuring uncertainties. However, looking at the distri- 
bution of the data points in the central 2.2 kpc of the galaxy, 
the new G emini data presented here, together with the Subaru 
data from Bre solin et al.|p010| ), display a much lower scatter 
(0.06 dex, 11 objects), compared to the RS08 data (0.21 dex, 
19 objects) pertaining to the same inner portion of the disk. A 
Kolgomorov-Smirnov test indicates that the probability of ob- 
taining the low scatter measured for the Gemini+Subaru sam- 
ple from the distribution of the RS08 data points is < 3%. 
This suggests that the dispersion in Fig. [3] might be related 
to the quality of the spectroscopic data drawn from the lit- 
erature, for example to the signal-to-noise (S/N) ratio of the 
[O III] A4363 line detections, rather than intrinsic azimuthal 
variations on scales of the order of a few hundred pc. It can 
also be seen that the RS08 abundances for the two innermost 
Hll regions (B43b and B 29) are significantly lower (~ 0.3 
dex) than those measured from the Gemini sample. 

To test the effect of the reported quality of the A4363 line 
detection on the abundance scatter, Hll regions for which 
this line was observed with a S/N < 6 were removed, together 
with 11 RS08 objects in common with the Gemini and Sub- 
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Table 2 

Reddening-corrected line fluxes and oxygen abundances 



ID 


[Oil] 


[Nem] 


[Oiii] 


Hi 


[Om] 


[Om] 


F(H/3) c(R/3) 


T e 


12 + log(0/H) 




3727 


3868 


4363 


4101 


4959 


5007 


(erg s _1 cm -2 ) (mag) 


(K) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) (9) 


(10) 


(11) 



9.. 
10. 
11. 
12. 
13. 
14. 
15. 
16. 
17. 
18. 
19. 
20. 
21. 
22. 
23. 
24. 
25. 



169 
213 
189 
245 
169 
195 
275 
181 
194 
192 
206 
239 
266 
197 
221 
178 
240 
238 
262 
242 
243 
219 
248 
224 
253 



± 17 
±20 
± 18 
±25 
± 16 
±20 
±26 
± 18 
±20 
±22 
±20 
±23 
±26 
± 19 
±21 
± 17 
±23 
±23 
±25 
±23 
±23 
±21 
±24 
±23 
±24 



4.7 ± 0.5 
7.2 ± 0.7 

4.2 ± 0.7 
4.9 ± 1.1 

8.3 ±0.8 
4.6 ± 0.8 
2.9 ± 0.7 

3.4 ± 0.4 



8.9 ± 0.9 
2.5 ± 0.3 
3.9 ± 0.7 
3.7 ± 0.6 
46.1 ±4.2 

6.4 ± 0.8 
27.7 ± 3.2 
49.4 ± 4.6 

8.1 ±0.8 

3.2 ± 0.4 

7.5 ± 1.0 



0.42 ± 
0.42 ± 



0.04 
0.07 



0.59 ± 0.08 



0.18 ±0.04 



0.48 ± 
0.30 ± 



0.12 
0.08 



2.78 



3.31 



±0.32 



±0.41 



23 ±2 
25 ±2 

25 ±2 

27 ±2 

28 ±2 

29 ±2 

26 ±2 
26 ±2 
26 ±2 
26 ±2 

25 ±2 

26 ±2 
26 ±2 
25 ±2 
25 ±2 

25 ±2 

26 ±2 

26 ±2 

27 ±2 
26 ±2 

25 ±2 

24 ±2 

26 ±2 

25 ±4 
25 ±2 



50 ±3 
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1.7 


xlO- 


14 


0.46 


7960 ± 300 


8.42 ± 0.09 


56 ±3 


171 ± 10 
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xlO- 


14 
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1.1 
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15 
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194 ± 11 
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0.67 
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17 ± 1 
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0.77 
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0.44 
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Note. — The line fluxes are in units of H/3 = 100. F(H/3) in column (8) is the measured H/3 flux, uncorrected for extinction. 
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Figure 2. Comparison between the data obtained for this work and results 
from the literature for: (a) the [O n] A3727 (circles) and [O in] A5007 (trian- 
gles) reddening-corrected fluxes (in units of H/3 = 100); (b) the [O in] A4363 
fluxes; (c) log(R.23); (d) 12 + log(0/H). The dashed line represents the one-to- 
one correlat ion line. The open and full symbols refer to the comparison with 
RS08 and Esteban et al. ( 2009] NGC 595 only), respectively. 

am observations done by the author. The result is shown in 
Fig. [4] The overall scatter is only slightly reduced to 0.13, 
still a factor of two larger than what is measured in the in- 
ner two kpc from the author's data alone. However, the ap- 



parent abundance peak shown by the RS08 data at 1.5 kpc 
from the galaxy center, as well as the lack of high O/H val- 
ues in the central kpc present in the RS08 sample, are now 
gone. The radial oxygen abundance gradient measured from 
this subsample (short-dashed line in Fig. [4} is: 



12 + log(0/H) = 8.48 (±0.04) - 0.042 (±0.010) R kpc (1) 

where Rk pc is the galactocentric distance in kpc. A 2 a itera- 
tive clipping was also applied to the data, providing a compa- 
rable result, but with smaller errors (long-dashed line): 



12 + log(0/H) = 8.50 (±0.02) - 0.045 (±0.006) R kpc . (2) 

Because of the removal of several low-metallicity H II regions 
in the central two kpc belonging to the RS08 sample the slope 
of the gradient given above is larger than the one determined 
by RS08 (-0.027 ±0.012 dex kpc" 1 ). The intercept value is 
also larger, 8.50 vs. 8.36. The measurement of the new Gem- 
ini T e -based abundances in the inner two kpc of M33 thus 
shows that ( i) it is possible that the real abundance scatter in 
this galaxy is not as large as claimed by RS08; (ii) the ra- 
dial abundance gradient in M33 behaves 'normally', as seen 
in other galaxies, with the metallicity increasing smoothly to- 
wards the center, wi thout indications of a central dip, which 
has been ascribed by Magri ni et al. (2010) to a selection bias. 
The interpretation given here is that it might have to do in- 
stead with errors in the measurements. Magri nTet al.| f2010) 
have also proposed that using only giant H II regions to de- 
termine the abundance gradient yields a steeper gradient than 
using the full sample, regardless of luminosity. This is not 
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Table 3 

Abundances for additional H II regions 



ID 


R 


12 + log(0/H) 




(kpc) 




(1) 


(2) 


(3) 


B72 


1.18 


8.46 ±0.09 


B90 


1.40 


8.47 ±0.08 


B302 


2.07 


8.44 ±0.08 



Note. — Data from the Subaru observations of 
|Bresolin etaTlpoTol . 

confirmed by the abundance analysis carried out in the cur- 
rent paper. From the same nine giant H II regions they used, 
the abundance gradient is: 



12 + log(0/H) = 8.50 (±0.10) - 0.040 (±0.024) R kpc , (3) 

which is essentially the same as found above. 

The latter result, based on only nine regions, is useful to of- 
fer an explanation as to how various authors have come up 
with different values for the slope of the gradient in M33. 
The 1 a range for the slope in Eq.|3] virtually encompasses all 



^ the central H II 
regions), several of which have been determined from a sim- 
ilarly small nu mber of ~10 Hll regions. (The central abun- 
dances from |Vilchez et al.fT9 88 are model-based, and lead to 
a steep gradient in the inner portion of the disk; however, the 
direct measurements based on the new auroral line detections 
do not support such high abundance values). A significantly 
larger number of data points, well distributed i n radius, is re- 
quired to reduce the unce rtainty in the slope ( [Dutil & Roy| 
|200T]|Bresolin et al.|2009) . 

As a final remark, the spatial distribution of the H II regions 
in the central parts of M33 is very similar between the Gem- 
ini and the RS08 samples, with the majority of the objects 
located along the spiral arms. Although it is implicit in the 
distribution of the Gemini and Subaru data points in Fig. H] 
it is worth emphasizing that pairs of H II regions with nearly 
coincident galactocentric distances, but situated almost dia- 
metrically opposite each other relative to the galaxy center 
have identical oxygen abundances, within the errors. This is 
the case for B 90 vs. B 50, B 72 vs. B 33b and B 302 vs. C 1 Ab 
(th e abundances for the three H II regions studied with Subaru 
by [Bresolin et al. 2010 are summarized in Table [3] for con- 
venience). This implies that the ISM of M33 is well mixed 
azimuthally over scales of at least 4 kpc. 

3.1. The abundance scatter determined from R23 

Despite the fact that the subsample of H II regions selected 
from the literature with a reported detection of [O III] A4363 
better than 6 a displays only a slight improvement in the rms 
scatter around the exponential gradient compared to the full 
sample, the much smaller dispersion from the Gemini+Subaru 
observations by the author suggests that the real abundance 
fluctuations in M33 are actually smaller and comparable to 
what is observed in other well-studied spiral galaxies. To test 
this assertion, and to show that it not simply a consequence 
of the relatively small number of objects considered (11) in 




Figure 3. Radial oxygen abundance gradient in M33, using [Oin]A 4363- 
based a bundances from this paper ( full circles), Bresolin et al. ( 2010] open 
circles), Rosolowsky & Simon ( 2008 squares) and other sources in the litera- 
ture (triangles). The least square tit to the full dataset is shown by the dashed 
line. Vertical segments connect objects in common between the Gemini and 
RS08 samples. 



years, from -0.012 ±0.011 dex kpc" 1 ( 


Crockett et al.||2006) + 


to -0.05 ±0.01 dex kpc -1 ( Vflchez et al. 


1988 after scaling to cm 




Figure 4. Same as Fig. [3] with sample reduced to the highest signal-to-noise 
[O in] A4363 detections, as explained in the text. The lines represent fits to 
all the data points shown (short dashes) and to a smaller sample obtained by 
means of a 2 a iterative clipping procedure (long dashes). 

the inner disk alone, the oxygen abundances for the H II re- 
gion sample shown in Fig. [5] were computed with the R23 
abundance indicator, which uses only strong emission lines, 
so that the resulting abundances are largely independent of 
the random errors affecting the much weaker [O III] A4363 au- 
roral lines, and which propagate through the chemical abun- 
dance measurement. Despite the large systematic uncertain- 
ties affecting R23 and other strong-line methods, as shown by 
the comparison of T e -based abunda nces with differen t calibra- 
tions of strong-line indicators (e.g. |Bresolin et al.|2009| ), the 
abundance scatter should be comparable to what high-quality 
auroral line detections indicate (somewhat larger dispersions 
can be expected from the dependence of R23 o n the ioniza- 
tion condition s of the different H II regions, see Kennicutt & 
Garnett 1996). In fact, as mentioned earlier, Bresolin et al. 
( 2009| measured the abundance gradient of NGC 300 from 
[O III] A4363 line detections in 28 H II regions, finding an rms 
scatter of 0.05 dex. Comparable values, between 0.06 and 
0.08 dex, are measured from the same data set when the abun- 
dances are obtained from a variety of strong-line abundance 
indicators. 



The R23 calibrati on of McGaugh (1991), as given in Kuzio 
de N aray et a"L] ( |2004| ), was adopted, since it has the advantage 
over other available calibrations of accounting for the nebular 
gas excitation, as given by the y = log(03/02) parameter, with 
03 = [O in] AA4959, 5007/H/3 and 02 = [0 11] A3727/H/3. The 
upper branch calibration was used, considering the published 
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[N II] /Ha values, and the position of the targets in a R23 vs 
O/H plot (presented in Sec. [4]). 

The O/H abundance ratio thus calculated is shown as 
a function of radial distance in Fig. [5] for the same H II 
regions presented in Fig. [3] A subsample, shown with 
full symbols, was createa by selecting objects with (a) 
S/N ([Om]A4363)>5; (b) [Nem] A3868/H/3 < 0.45; (c) 
[O ill] A5007/H/3 < 5; (d) logR 23 < 0.90. These criteria will 
be justified and discussed at length in Sec. [4] and serve to iso- 
late H II regions that are not affected by extremely hard ion- 
ization fields, or objects for which the R23 -based abundances 
could be very uncertain. H II regions that do not satisfy these 
conditions are shown either with open symbols (low S/N in 
the [O in] A4363 line), double open symbols (strong [Ne in] 
and [Olll] lines) and crossed symbols (large R23). 

The dashed line in Fig. [5] represents a least square fit to 
the data belonging to the subsample thus obtained, limiting 
the fit to galactocentric distances smaller than 5 kpc (to avoid 
possible biases when reaching down to metallicities near the 
turnover region of R23, as explained in Sec. [4]). The result is 
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Figure 5. Radial abundance gradient in M33 measured from the R23 indi- 
cator. Hn regions from the new Gemini observations (red circles), the Sub- 
aru data of Bresolin et al. ( 2010 bue circles), the sample of RS08 (squares) 
and other data irom the literature (triangles) are included. The dashed line 
represents the linear fit to the subsample defined in the text and shown by 
full symbols. Double-lined symbols are used for Hn regions with very high 
[Nem] and [Olll] fluxes. Crossed symbols are used for additional objects 
with log(R23) > 0-90. The remaining open symbols represent H II regions 
with low S/N in the [O in] A4363 line. 



12 + log(0/H) = 8.82 (±0.02) - 0.033 (±0.005) R kpc . 

The slope is slightly shallower than the value obtained from 
the auroral lines, although the two are still consistent within 
the 1 <j errors. It is worth recalling here the dependence of the 
slope of galactic radial abundance gradients on the method 
chosen to measure the oxygen a bundances for their H II re- 
gions (e.g. Bresolin et al. 2009). A s an example, applying 
the widely-used R2 3 calib rations by |Tremonti et al.| |2004 ) 
and |Zaritsky et"aL] ( |1994| ) yields slopes of -0.044 ± 0.006 
dexkpc -1 and -0.051 ±0.007 dex kpc -1 , respectively. More- 
over, as mentioned earlier, there are important systematic dis- 
crepancies in the metallicities derived from strong-line meth- 
ods calibrated from theo retical model gri ds (as in the case 
of the R23 calibration by|McGaugh 1991 adopted here) and 
the auroral-line method. This explains the large offset be- 
tween the intercepts of the regression lines shown in Fig. [3] 
and Fig.|5] However, this has virtually no effect on the deriva- 
tion of the abundance scatter. 

The most important result to note is that the rms scat- 
ter of the data shown with full symbols in Fig. B] is only 
0.05 dex, which is virtually the same value foundrrom the 
[O III] A43 63 -based abundances for the innermost H II regions 
observed with Gemini and Subaru. Changing the R23 calibra- 
tion has virtually no effect, as the rms scatter remains between 
0.05 and 0.07 dex throughout the disk o f the galaxy. Follow- 
ing RS 08, the least-square method by Akritas & Bershady 
(1996), which computes the intrinsic scatter of the data, was 
applied to the R23 -based abundances, confirming that the ob- 
served scatter is due only to measurement errors. The data 
therefore indicate the absence of real azimuthal variations in 
the oxygen abundances. 

It should be pointed out that the RS08 H II regions located 
in the central two kpc, which showed a very large scatter in 
the T e -based abundance plot of Fig. [3] are now in excellent 
agreement with the new Gemini data, and have a compara- 
bly small abundance scatter. It should also be added that the 
selection based on the S/N of the [O ill] A4363 line has lit- 
tle effect on the measured rms scatter, but it helps to retain the 
targets with the best determined values of the [0 11] A3727 and 
[O in] AA4959, 5007 lines, used to define the R 23 indicator. 

The main conclusion that can be drawn from the analysis 



of the R23 -based abundances, supported by the discussion of 
the T e -based results in the inner two kpc, is that the abun- 
dance scatter in M33 is small (0.05-0.06 dex), comparable to 
the measuring errors, and similar in magnitude to what is ob- 
served in other well-studied spiral galaxies, such as NGC 300. 
We can therefore speculate that the significant intrinsic scatter 
seen in the r e -based abundances of M33 from the RS08 sam- 
ple is due to measurement errors in the weak [O in] A4363 
auroral line. From Fig. [4] and [5] it is also possible to claim that 
the oxygen abundance reaches its peak at the center of the 
galaxy, as is typically the case for spiral galaxies, and there- 
fore that there is no evid ence for the off-cente red abundance 
distribution suggested by Magri ni et al.| ( |2010| ). 



4. SYSTEMATIC EFFECTS OF A HARD IONIZING RADIATION ON 
THE R 23 ABUNDANCES 

In Fig. [5] a number of targets were found to lie systemat- 
ically below the sequence of Hll regions defining the R23- 
based abundance gradient in M33. In particular, in the in- 
ner part of the disk the oxygen abundances of several H II re- 
gions are ~0.3 to 0.5 dex below the regression representing 
the abundance gradient, and appear to belong to a secondary 
sequence of data points that merges into the main one at a 
galactocentric distance of about 7 kpc. 

In order to interpret the origin for this discrepancy, and 
justify the exclusion of these targets from the R23 -based 
abundance gradient analysis, the T e -based oxygen abun- 
dances of the same Hll regions included in Fig. [5] are shown 
as a function of R23 in Fig. [6] In this diagram, a number 
of extragalactic Hll regions with high-quality metallicity 
determinations, based on the [Olll]A4363 auroral line, 
are also included (small dots). They were drawn from the 
following sources: 



NGC 300: 
M101: 
NGC 2403: 

IC 1613: 
Hn galaxies: 

M31, NGC 2363, NGC 1741 
NGC 4395, NGC 4861 



Bresolin et al. (2009) 
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Figure 6. Oxygen abundances derived from the [O in] A4363 auroral line as 
a function of log(R.23). A sample of H II regions with high-quality measure- 
ments in different spiral and irregular galaxies (small dots) have been used 
to schematically define the upper and lower branches of R23 (shown by the 
two intersecting curves). H II regions from the new Gemini observations (red 
circles), the Subaru data of Bresolin et al. (2010 bue circles), the sample of 
RS08 (squares) and other data from the literature (triangles) are included. 
The symbols used are the same as in Fig. [5] 

These data illustrate the well-known bi-valued nature of R23, 
consisting of two H II region sequences, intersecting around 
12 + log(0/H) = 8.0. The curves in Fig. [6] are polynomial fits 
shown to guide the eye in recognizing the two branches. The 
H II region data for M33 are plotted with the same symbols 
used in Fig. [5] It can be immediately seen that the major- 
ity of the most discrepant points in Fig. [5] identified by the 
large doubled-lined symbols, lie at large R23 values and above 
the intersection between the upper and lower branches in 
Fig. [6] Very large R23 values (log R23 > 1) are relatively rare 
among extragalactic H II regions, and are encountered in high- 
excitation nebulae with large ionization parameters. A frac- 
tion of such objects show spectral features associated with the 
presence of massive stars emitting partic ularly hard ion izing 
radiation, such as early WN (WNE) stars ( |Crowther|20Q7] ). It 
is thus not surprising that the H II regions where RS08 have 
identified nebular He II A4686 line emission (LP. = 54.4 eV) 
are located, together with additional discrepant points, in the 
upper right corner of the R23 vs. O/H plot. Invariably in this 
zone of the diagram, even when the He II emission remains 
undetected or the ionizing sources cannot be unequivocally 
identified, the high-excitation lines [Nelll]A3868 (LP. = 41 
eV) and [O in] A5007 (LP. = 35 eV) are particularly strong. 

4. 1. Hll regions with large [Ne III ] and [Olll] line fluxes 

The large double-lined symbols in Fig. [6] were used to iden- 
tify H II regions where Ne 38 68 = [Ne in] A3 868 /H/3 > 0.45 or 
O5007 = [O ill] A5007/H/3 > 5. The use of both line ratios is 
somewhat redundant, since the two are strongly correlated 
(the relative values of the criteria adopted for Ne3868 and O5007 
are consistent with the trend found between the two line ra- 
tios from the enlarged H II region sample analyzed here). The 
choice of the adopted selection criteria is purely empirical 
and somewhat arbitrary, and was guided by the knowledge of 
the line fluxes observed in nearby nebulae ionized by WNE 
stars. For instance, in RCW 5 a nd RCW48 (Milky Way ), and 
N57C and N138D,B (LMC), |Ken5cutt et al.| ( [2000] ) mea- 
sured Ne3868 > 0.56 and O5007 > 5.8, while, for example, for 
the giant H II region 30 Dor in the LMC, ionized by a cluster 
of O stars, Ne 3868 = 0.40 and O5007 = 4.3. 

Typical upper branch H II regions in other galaxies are ob- 
served with Ne3868 < 0-40 and O5007 < 4.0, and highly excited 
Hll regions, such as #15 (C lAb) and #19 (B 23) in Table[TJ 
easily stand out among objects with much lower Ne3868 and 



O5007 values. However, it should be recalled that high Ne3868 
and O5007 values, also exceeding the limits adopted above, can 
be encountered in lower branch H II regions (low metallici- 
ties), or near the turnaround between the two R23 branches. 
This explains why a few Hll regions at galactocentric dis- 
tances > 6 kpc were rejected in Fig.|5]by the criteria adopted 
above, even though their abundances can agree with the abun- 
dance gradient defined by the remaining regions. In fact, 
the Hll regions at large radial distances are near or at the 
turnaround, as Fig. [6] shows. For this reason the linear fit to 
the R23 -based abundances has been limited to galactocentric 
distances smaller than 5 kpc. 

It should be noted that not all of the discrepant objects in 
M33 contain known Wolf-Rayet stars, and also that the pres- 
ence of WNE stars in an extragalactic H II region does not 
necessarily imply peculiarly large Ne3868 and O5007 values. 
This occurs only if the nebular ionization is dominated by the 
radiation emitted by the WN star(s), and the dilution result- 
ing from the presence of O stars is small. This agrees with 
the observation that the total H ionizing flux output of the dis- 
crepant M33 Hll regions analyzed here, estimated from the 
Ha luminosity, is consistent with the ionization being pro- 
vided by only 1-2 WNE stars. What could be a little surpris- 
ing is finding such a large fraction (~ 20% of the nebulae 
with [O III] A4363 detections) of highly excited H II regions 
in M33. To the author's knowledge, for no other galaxy has 
an H II region population been insofar found with such a large 
percentage of nebulae having log(R23) > 1. On the other hand, 
the typical luminosity of H II regions selected in extragalactic 
abundance studies is generally significantly higher than that 
of the majority of those observed in M33. Consequently, their 
integrated spectra would be less affected by the presence of a 
few sources of hard radiation amid a large number of ionizing 
O stars. In fact, none of the M33 H II regions with a total lu- 
minosity L(HaQ > 10 38 erg s" 1 (adopting fluxes pu blished by 
Hodge et al.|[2502| and [Relano & Kennicutt||2009| ) for which 
a spectrum is available displays notable levels of excitation. 
The only (marginal) exception in t he upper branch reg ime is 
IC 131 (Ne 3 868 = 0.47, O5007 = 3.9, |Vflchez et al. 1988 ). The 
reason in this case is possibly connected to the fact that, be- 
sides containing three WN stars, this giant H II region has also 
a peculiarly hard and exte nded X-ray emission, whose source 
is at present unexplained (Tullm ann et al.|2 009). 

4.2. Excitation diagnostic diagrams 

To clarify the effect of the presence of WNE stars or 
other sources of hard ionizingradiation on the emission line 
spectra of Hll regions, Fig. n\ shows a [N II] A6583/H& vs . 
[O in] A5007/H/3 diagnostic diagram ( [Baldwin et"aL]|1981| ), 
in which M33 Hll region data from the literature (including 
the recent work on He Il-emitting H II regions by Kehri g"et al.| 
|2010| ) are shown together wi th a small number o f Galactic 
and ex tragalactic H II regions (Kennicutt et al. 2000 ; Bresolin 
|et al.||2009| ) containing identified WNE stars (star symbols). 
TEeWR ring nebulae G2.4+ 1.4 (W Ol), NGC 23 59 (WN4), 
NGC 3199 and S308 (WN5) from |Estebanet~aLl { [1993) are 
also shown (crosses). Additional Hll regio n data for spi- 
ral and irregular galaxies from the literature (Bres olin et al.| 

j2005 , and those shown as small dots in Fig. [6} are in- 
[ for comparison. Nebulae lying above the curve 1 



1999 
elude* 



comparison. JNebul ae lying above the curve (taken 
from Kauffma nn et al.|[2003| ) are considered to be ionized 
by a harder ionizing field than provided by O stars. In this 
plot the selection of objects with O5007 > 5 for the discrepant 
points in Fig. [6] corresponds to the shaded area, assuming that 
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Figure 7. The [Om]/H/3 vs. [Nn]/Hai diagnostic diagram for Hn regions 
in M33 from this work (triangles) and from the literature (full circles). Hn 
regions in the Mi lky Way and the L arge Magellanic Cloud ( Kennicut t et al.| 
2000), NGC 300 (Bresolin et al. 2009) and M33 (this work and Kehr ig et al | 
2010 1 hosting WN stars are included (star symbols; most are identified). WR 
ring nebulae from Esteban et al. (1993) are shown with crosses. A dditional 
nebulae contained in spiral and irregular galaxies are taken from Bresolin 
et al. (1999 2005) and other sources (open squares) to help discern the H II 
region excitation sequence. The curve is the boundary between star-forming 
and AGN-dominated galaxies defined by Kauffmann et al. ( 2003). The shad- 
owed area includes high-excitation H II regions in the upper R23 branch hav- 
ing [Om]/H/5 > 5. 

log([Nii] A65 83/Ha) > -1.2 on the upper branch ( [Kewley &] 
|Ellison|2Q08l ). 

It can be seen that the presence of WNE stars can shift H II 
regions into the grey area and above the boundary line. Still, 
there can be low-metallicity (small [N ll]/Ha ratio) nebulae 
containing WNE stars (e.g. MA 1 in M33) whose position in 
the diagram agrees with the excitation sequence defined by 
the control sample. Moreover, some of the high-excitation 
Hll regions have no clear detection of associated hot stars 



(e.g. H673 in M33, Kehri g et al-pOTO] ). For 13 objects in 
the Gemini sample the spectral coverage extended to suffi- 
ciently long wavelengths that the [N II] A5755 could be mea- 
sured. From the knowledge of the electron temperature or 
assuming T e = 10 4 K the strength of the [Nil] A6583 was es- 
timated using the temden program in IRAF. These objects 
were then added to Fig. [7] as tri angles. One of these , region 
B 18, contains a WN star (Massey & Johnson 1998 ) clearly 
detected in the Gemini spectrum from its broad HellA4686 
emission line. This H II region lies to the right of the bound- 
ary drawn in the diagram. For another, region B 23, situated 
into the shaded area, no possible m atch in the Wolf-Rayet cat- 
alog by Masse y & Johnson| ( 1998| ) could be found (the same is 
true for B 90, observed with the Subaru telescope), although 
its excitation, as measured from the [O lll]/H/3 line ratio, is 
higher than for B 18. 



The position of some of the deviating Hll regions can 
be linked to the presence of shocks, which can enhance the 
strength of the [N II] emission line. T his is almost certainly 
the case for region 11 in NGC 1232 (Bresolin et al. 2005, 
identified near the bottom of Fig. [7]), whose spectrum shows 
strong [O I] A6300 and [S 11] AA6717, 6731 emission (thus be- 
ing a supernova remnant candidate). For other Hll regions 
the situation is less clear. For example, H673 in M33 has 
moderately strong [N II] and [S II] fluxes, bu t no [O I] emis- 
sion has been reported by Kehri g~etaL] ( |20 1 0| ) . Nitrogen self- 
enrichment is also an option to explain the large [N II] flux, 
since these authors report log(N/0) = -0.74 ± 0.28, versus a 
typical log(N/0) = -1 .2 for H II regions in M33 (Bresolin et al. 
[2010] ). 

Fig. [7] shows that also Mun ch 1, an outlying object in 
M81, studied in some detail by |Garnett & Shields| ( [1987| ), is 
among the restricted number of high-excitation extragalactic 
H II regions contained within the shaded area (Ne3g68 = 0.46, 
O5007 = 6.9). Given its location in the BPT diagram, it is pos- 
sible that future high-sensitivity spectra of this region will 
reveal nebular He II and/or WN-related stellar emission fea- 
tures. 

4.3. Hll regions with large R23 values 

In Fig. [5] it was decided to remove from the final sample 
also the Hll regions for which log(R23) > 0.9. The crossed 
symbols in the plot represent the three H II regions that were 
removed based only on this condition (rather than poor S/N, 
or large [Nelll] and [Olll] line fluxes). The condition on 
R23 is added because the turnover in the oxygen abundance 
vs. R23 relation occurs around log(R23) — 1- It is well known 
that here the R23 indicator is a poor abundance diagnostic, be- 
cause a wide range of abu ndances corres ponds to a narrow 
range of indicator values (McGaugh 1991 ). Studies of ra- 
dial abundance gradients in spiral galaxies should therefore 
exclude H II regions in the turnover region in order to avoid 
abundance biases or large measuring uncertainties. It should 
also be noted that the selection based on R23 effectively re- 
moves virtually all of the targets with high [Ne III] and [O III] 
fluxes. The previous discussion, however, is helpful to under- 
stand what kind of objects systematically deviate from the R23 
vs. O/H trend of Fig. [6] 

In conclusion, Hll regions ionized by uncommonly hard 
ionizing radiation (such as that generated by the presence of 
early WN stars, but other poorly constrained processes are 
also at work) are shifted to higher R23 values at constant O/H 
abundance ratios compared to 'normal' H II regions in a di- 
agram like the one shown in Fig. [6j As a consequence, R23 
would systematically underestimate the abundances for these 
objects, as observed in Fig. [5] because the upper branch bends 
to smaller O/H values witn increasing R23 . This also ex- 
plains the secondary sequence seen at the bottom of this fig- 
ure, progressively diverging from the gradient regression line 
with decreasing radial distance (equivalent to increasing R23). 
The empirical emission-line criteria outlined above allow to 
discriminate quite effectively against such discrepant H II re- 
gions. 

5. SUMMARY 

New GMOS/Gemini spectra of 25 H II regions located in 
the central two kpc of M33 have been analyzed. The oxygen 
abundances have been derived from the [O in] A4363 auroral 
line for eight of the targets. The scatter in the radial abundance 
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gradient has been investigated by combining existing samples 
of H II regions in this galaxy with the new data. The main 
conclusions of this work are summarized as following: 

( 1 ) the scatter in the oxygen abundance in the central two kpc, 
as derived from the auroral line measurements for the new 
Gemini sample, is approximately 0.06 dex, much lower than 
the value of ~ 0.21 dex obtained in the sam e inner portion 
of the disk fro m the H II regions observed by |Rosolowsky &| 
|Simon| ( |2008] ). 

( 2 ) the oxygen abundances estimated from the R23 metallicity 
indicator for a large sample of H II regions assembled from 
the literature in combination with the new observations yield 
a comparably small scatter (0.05-0.07 dex) over the whole op- 
tical disk ofM33. This dispersion can be explained simply by 
the measuring uncertainties. 

(3) no evidence is found for significant intrinsic azimuthal 
variations in the Hll region metallicity distribution in M33 
on scales from ~100 pc to a few kpc, as well as for a dis- 
placement of the abundance peak from the galaxy center. This 
result rules out large anomalies in the mixing of the gas. 

(4) a considerable fraction of the M33 Hll regions with au- 
roral line detections show spectral features revealing sources 
of hard ionizing radiation (He II emission, large [Ne III] and 
[Oin] line fluxes). In some cases these can be identified 
with early WN stars, but several cases exist where no obvi- 
ous source has been identified. For these nebulae the oxy- 
gen abundances derived from R23 are systematically under- 
estimated by large factors. However, it is straightforward to 
identify these highly excited H II regions from their large R23 
values and their position in the excitation diagnostic diagrams. 
Thus, abundance gradients derived for other spiral galaxies 
hosting such objects would still be correct, once these pecu- 
liar nebulae are removed from the abundance analysis based 
on the strong-line metallicity indicator. 
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